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Abstract. I consider ageing behaviour in two exactly solvable reaction-diffusion systems. 
Ageing exponents and scaling functions are determined. I discuss in particular a case in which 
the equality of two critical exponents, known from systems with detailed balance, does not hold 
any more. Secondly it is shown that the form of the scaling functions can be understood by 
symmetry considerations. 



1. Introduction 

Ageing phenomena may occur in systems which are rapidly quenched into a region in parameter 
space with several competing stationary states. These phenomena have been studied quite 
extensively in systems with detailed balance such as simple magnetic systems, which are initially 
prepared in a high-temperature state and then quenched to or below the critical temperature 
Tc |Q|2|E3l2|l5IEi|- O ne typically considers the autocorrelation and the autoresponse function, 
for which one expects scaling behaviour in the ageing regime, that is for t,s and t — s large 
compared to some microscopic timescale: 

C(M) := (0(x,t)0(x,a)) ~ s- b f c (t/s) 

where </>(x, t) is the order parameter describing the system and h(x, s) is a small external 
perturbation, for instance a magnetic field, o and b are critical exponents and the scaling 
functions fc and Jr behave for large arguments as 



My) y ~°° y- Xc/z 



jRvy) ~ y 



(2) 



where A# and Ac are new exponents and z is the dynamical critical exponent. For systems with 
decorrelated initial conditions Xr = Ac has been found and the relation a = b has been confirmed 
at criticality. Notice that for systems with detailed balance the latter condition is necessary 
because the fluctuation-dissipation ratio has to hold when the system reaches equilibrium. 
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The approach of local scale invariance (LSI) had been proposed in IE] to understand the 
form of C(t,s) and R(t,s) on the basis of symmetry considerations. It turned out that for 
magnetic systems the form of R(t,s) can be completely fixed by symmetries whereas C(t, s) is 
fixed up to a scaling function. Monte Carlo simulations in the Ising model and the XY model 
showed excellent agreement with these predictions [01 110| I12j , whereas renormalisation group 
computations yielded a correction at two loops the origin of which is still not quite clear. 

In the systems considered so far detailed balance holds. However many realistic systems do 
not possess this property. It is therefore interesting to see what happens if this condition is 
relaxed. Typical systems which lack detailed balance are reaction-diffusion systems, where 
particles undergo diffusion on a lattice and in addition, there are particle creation and 
annihilation processes. Numerical work has been done on the (fermionic) contact process [14lll5j . 
showing that at criticality dynamical scaling holds for the response function and the connected 
correlator, but that opposed to the above mentioned magnetic systems a ^ b. Here, I shall ask 
the question: Are there exactly solvable systems without detailed balance where a ^ 6, and if 
so, can the form of the scaling functions still be understood with the help of the theory of local 
scale invariance? 

In this text I shall consider two specific models without detailed balance. On the one hand 
they allow for exact calculation of the ageing exponents and the scaling functions, on the other 
hand they can be described by a field-theoretical formalism and can therefore be attacked by 
the theory of local scale invariance. This paper is organised as follows: In section 2, I introduce 
the models and present the main results. In particular I show that in one case one indeed 
encounters a ^ b. In section 3, I look at the same models from the field-theoretical perspective 
and demonstrate that the form of the scaling functions can indeed be understood by considering 
the symmetries of the models. 

2. The bosonic contact and pair-contact processes 

2.1. The models 

I consider a d-dimensional cubic lattice with one sort of particles A. It is important that the 
system is bosonic, which means that there is no restriction on the number of particles on one 
lattice site. In what follows I shall consider two different models: 

• The bosonic contact process (BCPD): The particles undergo diffusion with diffusion 
constant D, that is, they may jump to a nearest-neighbor site with rate D. Furthermore 
single particles can disintegrate with rate A or produce offspring with rate /i. The newly 
produced particles then sit on the same lattice site as the original particle. 

A^A, A^%, A^2A 

This process has been used to model the clustering of biological organisms in • 

• The bosonic pair-contact process (BPCPD): Also in this process there is diffusion with 
constant D. In addition two particles on one lattice site can coagulate with rate A or 
give birth to a third particle with rate \i: 

A A, 2A A, 2A^U 3A 

It has been shown by Paessens and Schiitz |17j that this model can be solved analytically 
for any finite time. This is mainly due to a formal analogy to the spherical model, where 
the control parameter a defined below formally replaces the temperature. 



The state of the system is characterised by the number of particles on each site. I denote this 
by { n } = {• • • > n x> • • •}) where the non-negative integer ra x gives the number of particles on site 
x. The temporal evolution can be described by a master equation, which can be turned into 
a Schrodinger-type equation by standard techniques |18| I19| I2f)j . One then has annihilation 
and creation operators a(x) and a^(x) at each lattice site, the nonvanishing commutation 
relations of which are [a(x),a^(y)] = <5 x ,y Then two quantities are defined: the state vector 
\n) := IIx( a ^( x )) nx 1 0) , where |0) is the vacuum state representing the empty lattice, and the 
vector |-P(i)) := J2{ n } -P({ n }> 0l n )> where P({n},t) is the probability to find the system in the 
state {n} at time t. The latter quantity obeys the equation dt\P{t)) = — H\P(t)), where the 
Hamiltonian H is given in terms of annihilation and creation operators and can be found in ^7] 
for the cases at hand. Time-dependent observables are obtained by passing to the Heisenberg 
picture and the temporal evolution of an observable g{t) is then given by the Heisenberg equation 
of motion 

d t g(t) = [H,g(t)] (3) 

Finally, the average of g(t) is calculated as {g)(t) := (s\g(t)\P(0)), where (s\ is a coherent state 
vector with the properties (s|a^(x) = (s\ and (s\H = 0. The quantities of interest are: 

• The local particle density 

p(x, t) := (at(x, i)o(x, t)) = (a(x, t)) (4) 

where the special property of the state (s| has been used in the last equality. 

• The connected two-point correlator: 

G(x - y,t,s) := (a(x, t)a(y, s)) - (a(x, t))(a(x, s)) (5) 

Here and in what follows I assume spatial translation invariance, so that two-point quantities 
depend only on the difference of the spatial coordinates. If scaling behaviour is found, the 
corresponding scaling function will be denoted by fa in analogy to (0). 

• The response function: To compute this quantity, one adds a small perturbation 
J2 X M x ) t) a ( x ) to the Hamiltonian H. This corresponds to spontaneous particle creation 
at an empty lattice site with rate h(x,t). Then the response function is simply defined as 

\ 5(a(x, t)) , , 

= lifer (6) 

The Heisenberg equation of motion (j3J) is used to derive differential equations for the quantities 
For the particle density and the two-point correlator the following initial condition are 

chosen 2 

p(x,0)=p Oj G(x,0,0)=0 (7) 

whereas the response function is required to be a delta-peak at t = s. All these equations can 
be solved by standard techniques. 

2.2. Results 

For the particle density, one finds the following results for both processes I17j 

• For [i > A particle creation outweighs particle annihilation and p(x, t) diverges. 

2 This ensures that the system is translationally invariant. Furthermore it can be shown that this choice can be 
realised with an initial Poisson distribution of the particle density at each lattice site. 



• For fi < A particle annihilation is stronger and the systems runs into the empty lattice 
state. 

• Only if \x = A creation and annihilation of particles are of equal strength. In this case one 
has p(x, t) = po for all times. 

In the sequel I shall only look at the most interesting case when A = \i. Then the creation 
and annihilation processes balance each other, one says that one is on the critical line. As the 
particle density remains constant, one needs to look at the variance cr p (t) of p(x,t) for further 
insights. Notice that cr p (t) is equal to G(0, t, t) up to a constant. One finds the following results 

pmini 

• For the BCPD the correlator G(0, t, t) behaves as t~ 2 +1 (for t — > oo) and one has a diverging 
variance cr p (t) if d < 2 3 , otherwise the variance is bounded. This means that for d < 2 
diffusion can not spread particles evenly so that they accumulate on very few lattice sites - 
a clustering transition occurs 

• For the BPCPD the control parameter a is defined by 

a := fi/D. (8) 

It measures the strength of the creation and annihilation processes in comparison to the 
diffusion process. It turns out that there is a critical value ac > and the behaviour of 
G(0,t,t) depends on whether a is larger, equal or smaller than ac- More precisely 



a < ac, d > 2 


a = ac, d > 4 


a = ac, 2 < d < 4 


a > ac or d < 2 


G(0, t, t) < - z ^5 > const 


G(0,M) 


G(0,t,t) ^ti- 1 


G(0, t, t) 1 ~°° exp(t/r) 



0) 



For small a, diffusion is dominant and the system stays more or less homogeneous. But if 
a is large enough, there is again a clustering transition. 

Finally the two-time quantities [2^ are considered. For the response function scaling 
behaviour is found in the ageing regime with the result 

R(t,s) =r QS -i +1 ((t/s)-% +1 - l) (10) 

from which the critical exponents a and Ar can be derived, as well as the scaling function /jj(y) 
For the connected two-time correlator, one also finds the scaling behaviour © in the BCPD, 
and in the BPCPD for the first three cases discussed in Q. The ageing exponents can be found 
in table l[T|). whereas the scaling function is given by the integral expression 

f G (y) = g C de e^- h \ y + 1 - wyi (11) 

Jo 

Note that a ^ bis found at criticality 21 . This entails in particular that there is no non-trivial 
analogue to the fluctuation-dissipation ratio known from magnetic system and that it is not 
possible to define in a straightforward way an effective temperature characterizing the system 
as suggested in [2*2*] . 

3 In the case d — 2 there is a logarithmic divergence. 



3. Local scale invariance 



In this last section I consider the same processes from a field-theoretical point of view. First 
I define the fields </>(x, t) := a(x, t) — po and 4> := a^(x, t) — 1. In this way the correlator 
(4>(x, t)(j)(x', s)) equals the connected correlator ®. Then, by taking the continuum limit, the 
Hamiltonian H is turned into a field-theoretical action £[</>,</>] |2()| [23], from which n-point 
correlators can be computed via 

tn)exp(-H[<f>, </>]). (12) 

The action can be split up into two parts as £[</>, cj>] = En [4>, <$>} + T, no i se [(p, cj>], where the first part 
E [(f>,(j>] is given by 

E [(j),$\ = JdRdu [<f>(2Md u - V 2 )0] for the BCPD , . 

E [<M] = JdRdu [c^(2A^c\ -V 2 )0-a0 2 (/> 2 ] for the BPCPD. ^ ' 

Both fields under the integral depend on the integration variables R and u, but here and 





bosonic contact process 


bosonic pair-contact process 




a < ac 


a = ac 


a 


2 1 


2 1 


2 L 


b 


2 1 


2 L 


if 2 < d < 4 
| - 2 if d > 4 



Table 1. Ageing exponents of the critical bosonic contact and pair-contact processes in the 
different regimes. Xr = Xq = d and z = 2 was found for all cases. The results for the bosonic 
contact process hold for an arbitrary dimension d, but for the bosonic pair-contact process they 
only apply if d > 2, since ac = for d < 2. 

in what follows I shall suppress the arguments of the fields for space considerations. A4 is a 
parameter related to the diffusion constant. The exact form of the second part is somewhat more 
involved and can be found in [231125 • Here it suffices to state that £ n0 j se [<?!>, 4>\ has the form 
[<M= En. m <3 flnm / dHdu(f) n (f) m , where the constants a nm vanish whenever n > m. The 
point of this split-up is the following: It can be shown that in both cases, Xo[</>, 0] has nontrivial 
symmetry properties. This is a well-known fact for the BCPD 8^ as the corresponding evolution 
equation for eft is a free Schrodinger equation, but has only been shown recently [251 124j for the 
case of the BPCPD. From these symmetry properties the so-called Bargmann superselection rule 
can be inferred, stating that (cp ■ ■ ■ 4> 4> ■ ■ ■ 4>)q = unless n = m 4 . This entails that the response 

n m 

function is independent of S„ o j se [0, cp] [HI [21] and is given by 

R(x - x', t, s) = (0(x, t)0(x', s)) (14) 

where (• • -)o denotes the average with respect to Sq[0, <f>]- Further investigations reveal, that 
the connected correlator is essentially given by a sum of integrals over three- and four-point 

4 For the BCPD and an even number of fields, this can also be inferred from causality combined with Wick's 
theorem 



functions, which have also been computed with respect to £o[<^> 4 > ]- Therefore, the main task is 
to determine the two-, three-, and four-point functions of the theory described by So[0, <j>] only. 

Let us consider the BCPD first. In this case the evolution equation for the field <f> is a free 
Schrodinger equation. The symmetry group of this equation, i.e. the group of transformations 
carrying solutions to other solutions, is the well-known Schrodinger group 27 . An element g of 
this group acts on space-time coordinates as (x, t) — > (x', t') = g(x, t) with 

t^t=———, x^x = — — ; ad - fa = 1 (15) 

7* + o jt + o 

where 1Z ist a rotation matrix and a, /3, 7, S and v are parameters. A solution ^(x, t) of the free 
Schrodinger equation is transformed as 

¥(x,t) - CW(x,t) = fgig-i^Mg- 1 ^)] (16) 

where the companion function f g is known explicitely. A field transforming in this way is 
called quasiprimary and the important assumption is to identify the fields <j) and <j> as the 
appropriate quasiprimary fields of the theory. One can show that n-point functions build from 
quasiprimary fields satisfy certain linear partial differential equations involving the generators 
of the Schrodinger group. Solving these equations, one finds that the response function can be 
fixed completely with the result |241 I28| 



R ( X -x<,„)= ro(t - s )^> (lf'- n, ^(">Xf^l) (17) 



where x\ and xi are free parameters which can be adjusted so that this expression is in line 
with the result (| !()[). To deduce the exact functional form of R(x — x',t,s) is the actual pre- 
dictive power of this theory, providing new information which can for instance not be obtained 
by simple scaling arguments. Also three- and four-point functions can be fixed to a certain 
degree by the symmetries. There is, however, a degree of freedom that remains in the form 
of undetermined scaling functions. One can show that the result (|llj) can be reproduced by a 
suitable choice of these functions |24j . but I shall not give the results here for limitations of space. 

The procedure works similarly for the BPCPD. However, here the Schrodinger equation one 
has to consider is non-linear, which leads to a modification of the generators of the symmetry 
group, see |251 126j . The result for the response function in this case is |24| 

1+ s u ^(x,+x,) ft\^ (xi ~ X2) ( M(x-x') 2 \ T ft t-s a \ 
R(x - x', t, s) = r (t - (- ) exp j ) 

(18) 

with an arbitrary function ^ and another parameter y. Also this expression can be brought into 
agreement with (jlU|) . One can see that for the BCPD the symmetries fix completely the form 
of the response function whereas for the BPCPD an arbitrary scaling function remains. We can 
argue however, that the universal scaling function \& should not depend on the non-universal 
quantity a, which leads to the statement • ^rj^, (t- s )v ) = ^ OT another undetermined 

function This is in line with predictions obtained by renormalisation group arguments jSJ. 
Finally one can show that in the BPCPD also the result 1)11(1 can be reproduced correctly by 
appropriately adjusting the free parameters of the theory |24j . Again I can not treat the latter 
point more explicitly here but have to defer the reader to the quoted references. 



I conclude by summing up the main results of this paper. I have considered the bosonic 
contact and pair contact processes on the critical line. These are typical processes without 
detailed balance the absence of which leads to the violation of the relation a = b. However, also 
for these systems without detailed balance the approach of local scale invariance is still suitable 
to gain insight in the form of the scaling functions, even though the symmetries do not suffice, 
especially in the case of the pair-contact process, to fix completely the quantities of interest. 
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